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D I Abstract 



For a given barrier S and a one-dimensional jump-diffusion process X{t), starting from 

X < S, we study the probability distribution of the integral As{x) = /J'^^^^X(t) dt 
determined by X{t) till its first-crossing time ts{x) over S. In particular, we show 
that the Laplace transform and the moments of As{x) are solutions to certain partial 
• differential-difference equations with outer conditions. The distribution of the mini- 

mum of X{t) in [0,r5(2;)] is also studied. Thus, we extend the results of a previous 
paper by the author, concerning the area swept out by X{t) till its first-passage be- 
low zero. Some explicit examples are reported, regarding diffusions with and without 
jumps. 
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: 1 Introduction 
(N . 

CS| ' This paper deals with the first-crossing area, As{x) = /q^^^""^ X{t) dt, determined by a one- 
dimensional jump-diffusion process X{t) starting from x, till its first-crossing time, ts{x), 
over a threshold S > x, and it extends the results of a previous paper by the author ([2]), 
^ i concerning the area swept out by X(t) till its first-passage below zero. As for results about 
the integral of X{t) over a deterministic and fixed time interval, see e.g |5]. 
Notice that As{x) = X[o^rs{x)] ■ 'Ts{x), where X[o^rs{x)] denotes the time average of X{t) over 
the interval [0,r5(x)]. 

The first-crossing area has interesting applications in Biology, for instance in the frame- 
work of diffusion models for neural activity, if one identifies X{t) with the neuron voltage 
at time t, and ts with the instant at which the neuron fires, i.e. X{t) exceeds the potential 
threshold value S; then, As/ts represents the time average of the neural voltage till the 
first-crossing time over S. Another application can be found in Queueing Theory, if X{t) 
represents the length of a queue at time t, and one identifies the first-passage time ts over 
the threshold 5* with the overflow time, that is the instant at which the queue system first 
collapses; then, the first-crossing area represents the cumulative waiting time experienced 
by all the "users" till the congestion time. 
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As for an example from Economics, let us suppose that the variable t represents the quantity 
of a commodity that producers have available for sale and Y{t) describes the price of the 
commodity as a function of the quantity in a supply-and-demand model. Let us assume 
that Y{t) is bounded between and S > 0, and consider the distance of Y{t) from S, that 
is the process X{t) := S — Y{t). If we denote by r the amount of product at which Y{t) 
falls to zero (i.e. X{t) reaches the value S), then X{t)dt = tS — A, where A = Y{s)ds 
provides a measure of the total value that consumers receive from consuming the amount r 
of the product. 

In the present article, we complete the study carried out in [2]; in fact, for a one- 
dimensional jump-diffusion process, in place of the first-passage area below zero, we consider 
the analogous problem of first- crossing area over a positive barrier. 

Precisely, let be given a barrier S > and a one-dimensional jump- diffusion process X(t), 
starting from x < S; our aim is to study the probability distribution of the integral 
As{x) = X{t) dt determined by X{t) till its first -crossing time ts{x) over 5*. We 

improperly call As{x) "the first-crossing area of X{t) over S" . Indeed, the area of the 
plane region determined by the trajectory of X{t) and the t-axis in the first-crossing period 
[0,rs{x)] is J^^^""^ \X{t)\ dt which coincides with As{x) only if X{t) is non-negative in the 
entire interval [0,r(a;)]. 

Since the topic was studied quite extensively in [2] (though in the slight different situation 
of first passage below zero), we will omit some details. We will suppose that X(t) is the 
solution of a stochastic differential equation of the form: 



with assigned initial condition X{0) = x; here Bt is a standard Brownian motion, z/(-, ■) is 
a temporally homogeneous Poisson random measure (see Section 2 for the definitions), and 
the functions b{-),a{-),j{-, ■) satisfy suitable conditions for the existence and uniqueness of 
the solution (see Section 2). The coefficients regulate the drift (6), the diffusion (a), and the 
sizes of the jumps (7) which occur at (random) exponentially distributed time intervals. 
The process X{t) which is the solution of the equation (II. ip reduces to a simple diffusion 
(i.e. without jumps) if 7(x, m) = 0, and in particular to Brownian motion with drift /i if 
b{x) = /i and a{x) = 1. 



the first-crossing time over S of the process X{t) starting from x < S, and assume that ts{x) 
is finite with probability one. We will study the probability distributions of ts{x) and of the 
integral As{x) = J^^^^^ X{t) dt, as well as their moments; moreover, the distribution of the 
minimum of X{t) in [0,r5(a;)] will be studied. In particular, we will show that the Laplace 
transforms of As{x) and ts{x), their moments, as well as the probability distribution of the 
minimum of X{t), are solutions to certain partial differential-difference equations (PDDEs) 
with outer conditions. 

The paper is organized in the following way: Section 2 contains the statement of the 
problem and main results, in Section 3 some explicit examples are reported. Finally, Section 
4 is devoted to conclusions and final remarks. 




(1.1) 



Denote by 



mi{t > : X{t) > S\X{0) = x} 
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2 Notations, formulation of the problem and main re- 
sults 



Let X{t) G / := — oo <«</?< +00, a time-homogeneous, one-dimensional jump- 

diffusion process which satisfies the stochastic differential equation (SDE) (11.11) with assigned 
initial condition X{0) = x, where Bf is a standard Brownian motion and ■) is a Poisson 
random measure on (0, 00) x IR. Then, X{t) can be represented as 

i*t i*t /'+00 

X{t)=x + / h{X{s))ds+ I a{X{s))dBs+ / / -f{X{s),u)iy{ds, du) (2.1) 

Jo Jo Jo J~oo 

For the definitions of the integrals in the right hand side of (12. ip and the Poisson measure, 
see [12]. The coefficients 6,0" and 7 completely specify the law of X{t). In particular, an 
atom (t, u) of the Poisson random measure v causes a jump from x to x + 7(x, u) at time t, 
if X{t~) = X. We assume that u is homogeneous with respect to time translation, that is, 
its intensity measure E{v{dt^ du)) is of the form 

E[u{dt,du)] = dt7r{du) (2.2) 

for some positive measure vr defined on !B(1R) (S(1R) denotes the Borel a— field of subsets 
of IR), and we suppose that the jump intensity 

/+00 
7r{du) > (2.3) 
-00 

is finite. 

We make the following assumptions on the coefficients: 

Al 6, cr : J — )■ IR are continuous functions and a constant K > exists, such that, for every 
x,y e I : 

\b{x)-b{y)\<K\x-y\ 

b\x) + a'^{x) < A'(l +x2) 

A2 0" is a non-negative, bounded function and it is differentiable for every x belonging to 
the interior of /. Moreover, there exists a strictly increasing function G : IR"*" — )■ IR such that 
G(0) = 0, /q+ G-^{s)ds = +00 and 

\a{x)-a{y)\<G{\x-y\) 

for every x,y & I. 

Bl For every u and x E I : 

/+00 
7^(x,M)7r(c/M) < K{l+x^) 
-00 

B2 For every u and x,y E I : 




7(0;, u) — 7(2/, u)\'jT{du) < K\x — y 
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Remark 2.1 The conditions Al, A2, and Bl, B2 ensure that there exists a unique non- 
explosive solution of (11 .ip which is a temporally homogeneous Markov process (see [12], p!5]): 
A2 holds, for instance, if cr(-) is Lipschitz-continuous {G{x) = const -x), or Holder-continuous 
of order a G [1/2, 1) {G{x) = const ■ x"). 

Notice that, if 7 = 0, or = 0, then the SDE (II. ip becomes the usual Ito's stochastic 
differential equation for a simple-diffusion (i.e. without jumps). 

In the special case when the measure vr is concentrated over the set {mi,M2} = { — 1, 1} 
with TT^Ui) = 6i and '~f{ui) = ej, we can rewrite the SDE (II. ip as 

dX{t) = b{X{t))dt + a{X{t))dBt + esrfA^zl^) + eidNi{t) (2.4) 

where ei < 0, €2 > and Ni{t),t > are independent homogeneous Poisson processes of 
amplitude 1 and rates 61 and 62, respectively governing downward (A'^i) and upward (A^2) 
jumps. 

Let D be the class of function f(t,x) defined in IR"'' x I, differentiable with respect to t 
and twice differentiable with respect to x, for which the function f{t,x + 7(x,m)) — f(t,x) 
is TT— integrable for any {t,x). We recall the generalized Ito's formula for jump-diffusion 
processes giving the differential of a function f & D (see [T2]): 



df{t,X{t)) 



df df 1 d'^f 

-^{t,X{t)) + b{X{t))^{t,X{t)) + -a\X{t))^{t,X{t)) 

+00 



dt (2.5) 



+ ^(t, Xit))aiXit))dBt + I [fit, X{t) + 7(X(t), u)) - f{t, Xit))Hdt, du). 



—00 



The differential operator associated to the process X{t) which is the solution of (II. ip . is 
defined for any function f & D by: 

Lf{t,x) = L,f{t,x) + L^f{t,x) (2.6) 

where the "diffusion part" is 

1 d'^f df 
Ldf{t,x) = 2^'^(^)q^(^^^) + ^(^)^(^'^) 

and the "jump part" is 

r+00 

Ljf{t,x)= / [f{t,x + 'y{x,u)) - f{t,x)]7r{du). 



* 'df 



L 



^{s,X{s)) + Lf{s,X{s)) 



ds . (2.7) 



Then, from (12. 5p . taking expectation, one obtains 

E[f{t,Xm = fiO,X{0)) + E 

For a barrier S and x < S*, we define: 

Ts{x) = mi{t > : X{t) > S\X{0) = x} (2.8) 

that is the first-crossing time over S of X^{t) (i.e. the process X{t) starting from x) and 
suppose that ts{x) is finite with probability one. Really, it is possible to show (see [2S], [S]) 
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that the probabihty Pq{x) that X^(t) ever leaves the interval (— oo, S*) satisfies the partial 
differential- difference equation (PDDE): 



Lpo = 



(2.9) 



with outer condition: 



Po{x) 



liix> S. 



The equality po{x) = 1 is equivalent to say that ts{x) is finite with probability one. For 
diffusion processes without jumps (i.e. 7 = 0) sufficient conditions are also available which 
ensure that Ts{x) is finite w.p. 1, and they concern the convergence of certain integral 
associated to the coefficients of (11.11) (see Section 3.1 and also [12] , Hi 



Remark 2.2 Let us consider the special case when X{t) = X(t) + J(t), where X is a 
simple-diffusion (i.e. without jumps) and J is a pure-jump process, set m{t) = E{X(t)), and 
suppose that 

3 t > : m{t) + E{J{t)) > S, for any t > i; (2.10) 

then, P{ts{x) < oo) = 1. Indeed, P{ts{x) = oo) > implies that, for a set of trajectories 
having positive probability, it results X(t) < S for any t > 0. Thus, by taking expectation 
one obtains E{X{t)) = rri{t) + E{J{t)) < S for any t > 0, which contradicts fl^HID . 
If J{t) = eiNi{t) + e2N2{t) (cf. ([22D), condition fl2:T0D becomes m(t) + (ei^i + £3^2)^ > S; 
in particular, for €2 = e = — ei it writes m(t) > S + e{6i — 6*2)^. 

Let f/ be a functional of the process X; assume that ts{x) is finite with probability one, 
and for A > denote by 



MuAx) = E 



-A/;sWt/(x(.))ds 



(2.11) 



the Laplace transform of the integral J^'^^^^ U{X{s))ds. Then, the following theorem holds 
(we omit the proof, since it is quite analogous to that of Theorem 2.3 in [2]): 

Theorem 2.3 LetX{t) he the solution of the SDE (11. ip . starting from X [Q) = x < S; then, 
under the above assumptions, Mu^\{x) is the solution of the problem with outer conditions: 



LMu,x{x) = XU{x)Mu,xix), x <S 
Mu,x{y) = 1, ioTy>S 
lim^^_oo Mu,x{x) = 

where L is the generator of X, which is defined by ( 12. 6p . 



We recall that the n-th order moment of J^^^^^ U{X{s))ds, if it exists finite, is given by 
(n = l,2,...) : 



(2.12) 



□ 



Tn{x) 



E 



U{X{s))ds 











= (-1)" 









A=0 
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Then, taking the n-th derivative with respect to A in both members of the equation fl2.12p . 
and calculating it for A = 0, one easily obtains that the n-th order moment T„(a;) (n = 
1, 2, ... ) of Jq^^^^ U{X{s))ds, whenever it exists finite, is the solution of the PDDE: 

LT„(x) = -nf/(x)T„_i(x), x<S (2.13) 

which satisfies 

T„(x) =0, for X > ^ (2.14) 
and an appropriate additional condition. 

Indeed, the only condition T„(x) = for a: > S* is not sufficient to determinate uniquely the 
desired solution of the PDDE fl2.13p . because it is a second order equation. We will return to 
this problem when we will consider some explicit examples. Note that for a diffusion without 
jumps (7 = 0) and for U{x) = 1, f l2.13p is nothing but the celebrated Darling and Siegert's 
equation ([H]) for the moments of the first-passage time, and f l2.14p becomes simply the 
boundary condition T„(5') = 0. 

Remark 2.4 In certain cases, the first-crossing time and the first-crossing area oi X{t) over 
S can be expressed in terms of the first-passage time and the first-passage area of a suitable 
process below zero. For instance, let X{t) = x + Bf + fit, with /i > 0, i.e. Brownian motion 
with positive drift /x; then ts{x) = inf{t > : X{t) > S} has the same distribution as 
inf{t > : S — X + Bf — fit < 0} = t{S — x), where T{y) denotes the first-passage time of 
the process Y{t) = y + Bt — fit below zero. Moreover, it is easy to see that 

As{x) = [ X{t)dt = S-r{S-x)- A{S - x) (2.15) 
Jo 

where A{y) denotes the area swept out by Y{t) till its first-passage below zero. 
Thus, the moments of ts{x) can be soon obtained by those of t{S — x), while by using fl2.15p . 
the mean of As{x) can be obtained in terms of E(t{S — x)) and E{A{S — x)) (the first two 
moments of r and A were obtained in [2] ). Notice that, the computation of E{As{x)'^) 
requires instead also the knowledge of the covariance of r and A. 

Distribution of the minimum of X(t). 

Now, we will study the probability distribution of the maximum downward displacement 
(i.e. the minimum) of the jump-diffusion X{t) starting from x < S, till its first passage over 
S, that is, m{x) = m.mt(z[Q^Tsix)]{X {t)\X (0) = x}. For any z < x < S the event {m(x) > z} 
is nothing but the event "X(t) first exit the interval [z, S] through the right end S"; so, by 
the well-known result about the exit probability of a jump-diffusion from the right end of an 
interval (see [S]), we obtain that v{x) := P{xn{x) > z), as a function of x, is solution of the 
equation Lv = with conditions v{y) = 1, y > S; v{y) = 0, y < z. Thus, we get: 

Proposition 2.5 P(ja{x) < z) is the solution of the problem with outer conditions: 

Lw{x) = 0, X G {z, S) 

w{y) = 0,y>S (2.16) 

wiy) = 1, y < z 

□ 
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3 A few examples 



In this section we will compute explicitly the moments of ts{x) and those of the first-crossing 
area As{x) for certain jump-diffusion processes. We start with considering diffusions without 
jumps. 



3.1 Simple diffusions (i.e with no jump) 

Let X[t) be the solution of (11.11) . with 7 = 0, that is: 

idX{t) = b{X{t))dt + a{X{t))dBt 



\X{0) 



X 



(3.1) 



In this case ts{x) = inf{t > : X(t) = S\X{0) = x}, that is the first-passage time of X^{t) 
through S. 

Let us consider the functions (c is a constant) : 

0(x)=exp(^-^ ^ds^ (3.2) 

r 2 

^(x) = 0(x) / ds. (3.3) 

As it is well-known (see e.g [I2]), a sufficient condition in order that ts{x) is finite with 
probability one, namely the boundary S is attainable, is that the function ^(x) is integrable 
in a neighbor of S. 

Since the generator L coincides with its diffusion part Ld, by Theorem l2.3l we obtain that, 
for X < S, Mu^xi^) = E e~'^-^o^ u{x(s))ds jg solution of the problem with boundary 
conditions (M' and M" denote first and second derivative with respect to x) : 



x)M^,(x) + bix)Mi,ix) = \U{x)Mu,x{x) 



Mu,x{S) = 1 . (3.4) 

\im^^_^ Mu,x{x) = 

Moreover, by (I2.13p . (I2.14p the n-th order moments T„(x) of J^^^^ U{X{s))ds, if they exist, 
satisfy the recursive ODEs: 

-a\x)T;^{x) + b{x)T^{x) = -nU{x)Tn-i{x), ior x < S (3.5) 
2 

with the condition Tn(S) = 0, plus an appropriate additional condition. 

Finally, as regards the minimum m(x), it turns out that its distribution F{z) = P(xn{x) < z) 

is the solution of the problem with boundary conditions: 

Lw{x) = 0, X G {z, S) 

w{z) = 1 . (3.6) 

w{S) = 
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Example 1 (Brownian motion with drift /i > 0) 

Let he X(t) = x + fit + crBf, with a > and < a; < 5*. Without loss of generality, we 
can assume a = 1 (otherwise, dividing by a one reduces to this case). Note that, since the 
drift is positive, ts{x) = Tg{x) = inf{t > 0\X{t) = S} is finite with probability one, for any 

(x) 

X < S. Taking b{x) = a{x) = 1, the equation in ([MD for Mu,xix) = E{e~^fo^ ^ uix{s)ds-^ 
becomes 

^M^_,(x) + /xM^,,(x) - \U{x)MuAx) = 0. (3.7) 

(i) The moment generating function of Tg{x) 

By solving (13. 7p with U{x) = 1, we explicitly obtain: 

Mu,x{x) = cie^'i" + cse^^" 

where pi = — /i — -^//^ + 2A < 0, P2 = — + v^/i^ + 2A > 0; the constants Ci and C2 must 
be determined by the boundary conditions. Indeed, M[7a(— oo) = gives Ci = 0, while 

Mu,xiS) = 1 implies Ca = e^^^"^^^^)^. Thus, we get: 



Mu,x{x) = E{exp{-Xr^{x)) = exp[(/i - ^/J^+2X){S - x)]. (3.8) 

This Laplace transform can be explicitly inverted (see [H]), so obtaining the well-known 
expression of the density of Tg{x) : 



For /i > the moments T„(x) = i?(r^(x))" of any order n, are finite and they can be easily 
obtained by calculating { — l)"-[d"'Mu,x{x)/dX^]x=o- We obtain, for instance: 

E{r^,{x)) = E{{r^si^)r) = ^ + ^-^^ (3.10) 

(cf. Remark 12.41 and the results for BM with negative drift in p] ). Note that E{Tg{x)) — 
+00, as /i — > 0. As easily seen, for any x < S* it results Mu^x{x) — 1, as /i — )■ +oo, or, 
equivalently, t^(x) converges to in distribution, and so _E'((r^(x))") — 0, as /i -> +oo, for 
any n and x < S. 

(ii) The moments of A^g{x) = J^^ ^ {x + pt + Bt)dt 
For U{x) = X the equation (13. 7p becomes 

^Mi}^, + fiMlj^, = \xMu,x (3.11) 

with conditions Mu,xiS) = 1, Mu,x{-oo) = 0; now Mu,x{x) = E(e"^^s(^)). Unfortunately, 
its explicit solution cannot be found in terms of elementary functions, but it can be written 
in terms of the Airy function (see [18], [13]) though it is impossible to invert the Laplace 
transform Mu^x to obtain the probability density of Ag{x). 

In the special case /i = 0, it can be shown (see [TS], [2]) that the solution of (13. lip is: 

Mu,x{x) = 3'/'r ('^^ Ai(2i/3A^/='(5 - x)) (3.12) 



where Ai{z) = vr -^^y z/3 -ft'i/s (f-^^^^^) denotes the Airy function, and K^{z) is a modified 
Bessel function (see [1]). Calculating the derivative with respect to A in fl3.12p . we obtain: 



By using the fact that Ai'(O) = - (S^/^F (|))"^ (see P), it follows that -§^Mu,x{x)\x=o = 
— oo, namely, the expectation of ^^(x) is infinite. Notice that, unlike the Laplace transform 
inversion reported in [2] for a similar case (see equation (3.12) therein), we are not able 
to invert explicitly the Laplace transform f l3.12p to find the density of A^; in fact, in the 
case considered in [2] this was possible, thanks to an integral identity (see e.g. [IB], [13] ) 
involving the modified Bessel function, which serves the purpose when the support of the 
candidate density of is (0, +oo). 

For yU > 0, we will find closed form expression for the first two moments of A'^{x), by 
solving f l3.5p with U{x) = x. For n = 1, we get that Ti{x) = E{Ag{x)) must satisfy the 
equation: 



\Tl\x) + fiTlix) = -X 
Ti{S) = 



, ^ ^^^^^^ 



Of course, the only condition Ti[S) = is not sufficient to uniquely determinate the 
solution. The general solution of fl3.13p involves arbitrary constants Ci and C2 and, as easily 
seen, it is given by 



2 

Ti(x) = ci + C2e-2^^-^+ 



2/i ■ 2^2 • 

By imposing that Ti{S) = and that for any a; < it must be Ti{x) — > 0, as /i — +00, we 
find that the mean first-crossing area is 



S + X-- 



(3.14) 



2/i 

This formula also follows by taking expectation in f l2.15p . and using that Eijiy)) = -, 

E{A{y)) = f^ + ^{seem)- 

As far as the second moment of Ag{x) is concerned, we have to solve (13. 5p with U{x) = x 
and n = 2, obtaining the equation for T2{x) = E[{Ag{x)Y] : 



ln\x) + fiTiix) = -f (5 -x){S + x- l//i) 
T2{S) =0 



As before, the only condition T2{S) = is not sufficient to uniquely determinate the solution. 
The general solution of (13.150 . which involves two arbitrary constants Ci and C2, is given by 

T2{x) = ci + C2e-2'^^ + Ax^ + Bx^ + Cx^ + Dx, 

where 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Figure 1: Estimated density of the first-passage area Ag{x) oi X(t) = x + fit + Bt over S, 
for S* = 2, X = 1 and several values of yU. From top to bottom, with respect to the peak of 
the curve: /i = 3; = 2; = 1.5; /i = 1.2; = 1. 

By imposing that, for any x < it must be T2(x) — )■ 0, as /i — > +oo, we find C2 = 0; 
moreover, by T2{S) = we get ci = —{AS'^ + BS^ + CS*^ + DS). Thus the second order 
moment of the first-crossing area is 

E[{A^s{^)f] = ^i.^^ - S^) + B{x^ - S^) + C{x^ - S^) + D{x - S), (3.16) 

where the constants A,B,C,D are as above. Finally, by (13.141) and (I3.16p . one can obtain 
the variance of Ag(x). 

Notice that £'[(yl^(x))^] has to be the only non-negative solution to (13.151) . Indeed, expression 
(I3.16P loses meaning if it becomes negative for some choice of x, fi and S; in that case, the 
second moment of Ag{x) does not exist. 

Since a closed form expression for the density of the first-crossing area Ag{x) cannot be 
found for /i > 0, it must be obtained numerically. As in the case of the first-passage area 
below zero (see [2], [3]), we have estimated it by simulating a large number of trajectories of 
Brownian motion with drift /i > 0, starting from the initial state x > 0. The first and second 
order moments of the first-crossing time Tg and of the first-crossing area Ag thus obtained, 
well agree with the exact values. For S* = 2, x = 1 and several values of /i > 0, we report in 
the Figure 1 the estimated density of the first-crossing area, regarding X{t) = x + fit + Bt. 
For some values of parameters we have compared the estimated density of the first-crossing 
area with a suitable Gamma density; in the Figure 2, we report for 5 = 2, x = 1 and 
/i = 1.5, the comparison of the Laplace transform Mx{x) = E[e~^'^si^)'j of the estimated 
first-crossing area density and the Laplace transform of the Gamma density with the same 
mean and variance. Although the two curves agree very well for small values of A > 
(this implying a good agreement between the moments), for large values of A the graph of 
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Figure 2: Comparison of the Laplace transform Mx{x) = E{e"^^s(^)) of the estimated density 
of the first-passage area Ag{x) (lower curve) and the Laplace transform of the Gamma density 
with the same mean and variance (upper curve), as functions of A, for 5 = 2, x = 1 and 
/i = 1.5. 



the Laplace transform of the estimated density of Ag{x) lyes below the other one, which is 
compatible with an algebraic tail for the distribution of Ag{x). 

(iii) The distribution of the minimum m^(x). 

Set 



m^(x) 



min \x + lit + BA. 



Its distribution function P{m^{x) < z) = w{x) is the solution of the problem with boundary 
conditions: 



^w"{x) + fiw'{x) = 
w{z) = 1, w{S) = 

By solving the above equation, we obtain, for z < x : 



P{m''{x) <z) = w{x) 



(3.17) 



(3.18) 



g— g— 2/iS 

Then, calculating the derivative with respect to z, we get the probability density of m^(a;) 



2/16-2'^^ (( 



-2iiS 



(—00,2;] (2^) 



(3.19) 



For /i = 0, taking the limit in the above expression, we get /niO(z)(^) = (s-z^ l(-oo,a;](2;), and 
so the moments of xn^(x) of all orders n > 1 are infinite. 
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On the contrary, for /i > the minimum m^(x) possesses finite moments of all orders, 
because fmf^{x){z) behaves like const / e~'^'^'' , for z — — oo. 

Definition 3.1 We say that a one- dimensional diffusion X{t) with X{0) = x, is conjugated 
to BM if there exists an increasing differentiable function u{x), with u{0) = 0, such that 
X{t) = u-^ {Bt + u{x)) . 

Remark 3.2 If X{t) is conjugated to Brownian motion via the function u, then for x < S : 
Ts{x) = mi{t > : X{t) = S\X{0) = x} = M{t > : u{X{t)) = u{S)\u{X{0)) = u{x)} 

= inf{t >0:Bt + u{x) = u{S)} = r|(M(x)) 

where S' = u{S) and Tg,{y) is the first hitting time to S' of y + Bt (i.e. BM without drift 
starting from y). Thus, the first-passage time through 5" of the process X(t) starting from 
X < S", is nothing but the first hitting time to S' of BM starting from u{x), and 

rrsix) r-r^iiuix)) 

As{x)= / X(t)dt= I u~^{Bt + u{x))dt. 

Jo Jo 

Notice that, though ts{x) turns out to be finite with probability one, it results E{ts{x)) = 
+00. Moreover 



Mu,x{x) =E exp ( -A / U{X{t))dt \ = E 
where V{s) = U (n~-^(s + u{x))) . 

Therefore, the Laplace transform Mu^x{x) of Jq'^''^' U{X{t))dt, associated to the functional U 

of the process X, is nothing but the Laplace transform My^{y) of Jq^'^^'' V{Bt)dt, associated 
to the functional V of BM, where y = u{x). 

This means that the equation (13.41) is easily reduced to the analogous equation for BM 
starting from u{x), with U replaced by V. 

Example 2. 

A class of diffusions conjugated to BM is given by processes X(t) which are solutions of 
SDEs such as 

dX{t) = ^a{X{t))a'{X{t))dt + a{X{t))dBt, X{0) = Xq (3.20) 

with (t(-) > 0. Indeed, if the integral v{x) = f^^ '^r)^'^ convergent for every x, by Ito's 
formula, we obtain that X{t) = v~^{Bt + v{xq)). 

(i) (Feller process or Cox-IngersoU-Ross (CIR) model) 
For a, 6 > 0, let X{t) be the solution of the SDK 



T^,(U(X)) 

exp I -A / V{Bt)dt 



r-Ts{x) 



dX{t) = {a + hX{t))dt + ^/X{t) V dBt , X(0) = Xq > (3.21) 

(note that, although y/x is not Lipschitz-continuous, the solution is unique because ^/x is 
Holder-continuous of order | (see condition A2)). The process X{t) turns out to be non- 
negative for alH > (see [7], [9]). If 6 = and a = i, the SDK (I3A) becomes: 



dX{t) = -dt + ^/X{t) VO dBt , X{0) = xq 
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and X{t) turns out to be conjugate to BM via the function v{x) = l^fx i.e. X(t) 
\{Bt + 2^)2; the SDE is obtained by taking a{x) = VxVO in (KW^ . 

(ii) (Wright & Fisher-hke process). 
The diffusion described by the SDE: 



dX{t) = {a + hX{t))dt + ^JX{t){l- X{t))yQ dBt, X(0) = xq G [0, 1] 



(3.22) 



with a > and a + b < does not exit from the interval [0, 1] for any time (see [7], [9]). 
This equation is used for instance in the Wright-Fisher model for population genetics and in 
certain diffusion models for neural activity [22] • For a = \ and b = X(t) turns out to be 
conjugated to BM via the function v{x) = 2arcsinA/x, i.e. X{t) = svo?[Bt/2 + arcsin^/xo). 
For these special values of parameters, the SDE fl3.22p becomes: 



dX{t) 



\ - \xit) ) dt 



'X{t){l- X{t)) dBt ,X(0) = xo e [0,1] 



and it is obtained from fl3.2Up by taking o"(x) = 
^ < 1, it results As{x) = X{t)dt < rs(x), since < X(t) < 1. 

Example 3 (Ornstein-Uhlenbeck process) 
Let X{t) be the solution of the SDE: 



xd - x) V . Notice that, for < x < 



dX{t) = -fxX{t)dt + adBt, X(0) 



X 



(3.23) 



where /i and a are positive constants. By calculating the functions and ^ in (13. 2p and 
(13. 3p . we obtain: 



and 



e(x) 



0(x) = const ■ exp 



2 RF ( ux^ 
const ■ — w — exp — - 



xv2/i 



cr 



where $ denotes the distribution function of the standard Gaussian variable. Since .^(x) is 
integrable in a neighbor of x = 5, the boundary 5* is attainable. 

The exphcit solution of (l3:23|) is X{t) = e"/^* (x + J^^aef^'dBs^ (see e.g. g], [6]). By using 

a time-change, we can write /J ae^'^dBg = Bp{t), where pit) = ^ [e^^* — 1) . Then, one gets 
X{t) = e~'^*(x + B{p{t)). Let us consider e.g. the time dependent boundary S(t) = H(t) = 
ae~'^*, then: 

Th{x) = inf{t > : X{t) = H{t)\X{0) = x} = inf{t > : B{p{t)) + x = a} 



and so 



that is Th(x) 



p{th{x)) = mf{s > : X + Bs = a} = r^(x 



P~\t^{x 



;)), where r^(x) is the first-passage time through a of Brownian 



motion starting from x. Thus, th{x) has density 

frHix){t) = frB(^,^){p{t))p'{t) 
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where /rB(a;) denotes the density of r^(x), which is given by (13.91) with fi = and S = a. 
Notice that even the mean of th{x) is impractical to be directly calculated by using fTH(x)if)- 
However, for a constant boundary S it results (see e.g. [21]): 

E{ts{x)) = i{v/7r [01 {S^/Ji/a) - 0i {x^/Jl/a)] + (S^/JI/a) - ipi (x^/JI/a) }, 



where 



0lW 



°° ^2fc+l 



fc=0 



{2k + l)k\' 



Notice that, if S* > 0, as seen by using Hospital's rule, E{ts{x)) tends to infinity, for /i — )• 0, 
since the OU process reduces to BM; also for fi — > +oo it tends to infinity, which is indeed 
natural, since the drift tends to — oo. 

As far as the minimum m(x) is concerned, by Proposition 12.51 its distribution w{x) = 
P(m(x) < z) satisfies the problem: 

|(T^w"(x) — fixw'{x) = 0, 2; < X < 5* 



w{z) = 1, w{S) = 
whose solution is: 

P{m{x) <z) = (^j^ e^'^dt 
Thus, the density of m(x) is: 



e^* dt 



z < X. 



d 
dz 



P(m(x)) < z) 



e5^' J e^^^dt e^^^d?j 



L(-oo 



3.2 Diffusions with jumps 



Example 4 (Poisson process) 

For X > 0, let us consider the jump-process X{t) = x + Nf, where Nt is a homogeneous 
Poisson process with intensity 9 > 0, namely Nq = and its jumps, of amplitude 1, occur 
at independent instants, exponentially distributed with parameter 6. This means that 

PiNt = k) = forfc = 0,l,... 

The infinitesimal generator of the process is defined by: 

Lg{x) = 9[g{x+l)-g{x)], geC'{JR) 
and, for 5* > and x < S, Ts{x) = mi{t > : x + Nt> S}. 
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By Theorem 12.31 with U{x) = 1, it follows that the Laplace transform Mu^x{x) of ts{x) 
is the solution of the equation LMu^x{x) = AM{/a(x), with outer condition Mu^x{y) = 1 for 
y > S. By solving this equation, we get: 

. ^ f(«Tl)^~" if^-xGlN 
Xiwx) if 5 -x^ IN 

where [a] denotes the integer part of a. Notice that the condition M[/a(— C)o) = also 
holds. Thus, recalling the expression of the Laplace transform of the Gamma density, we 
find that ts{x) has Gamma distribution with parameters {S — x,6) if 5" — x is a positive 
integer, while it has Gamma distribution with parameters {[S — x] + 1,6) if S* — a: is not an 
integer. This fact also follows directly, if one considers the nature of the Poisson process. 
The moments Tn{x) = E [Tg{x)] are soon obtained from the density or also by the formula 
T„(x) = (-1)" [§,MuAx)\^^^.V^e have: 



E{ts{x)) 

and 




if 5 - X G IN 
if 5 - X ^ IN 



(^ + ^ if^-xelN 
+ if5-x^lN 



^ if 5 - X G IN 



Thus: 

Varirsix)) = { 

K JJ ) [s~^ if 5-x ^ IN 

By Theorem 12.31 with U{x) = x, we get the Laplace transform Mu^\{x) of As{x) as the 
solution of the equation LMu^\{x) = XxMux^x), with outer condition Mu^x{y) = lioi y > S. 
By solving this equation, we get: 

M ( \ - l^^"" ' + + + 1)) ■ • • + A(x + {S-x)- l))}-i if 5 - X G IN 

uA^) - . ^ xx){e + A(x + 1)) ■ ■ ■ + A(x +[S- x]))}-i if 5 - X ^ IN 

Notice that the condition Mu^—oo) = is fulfilled. 

We observe that Mjjxix) turns out to be the Laplace transform of a linear combination 
of S* — X independent exponential random variables with parameter 6, with coefficients x, x + 
1, . . . , — 1, if S* — X is an integer, while it is the Laplace transform of a linear combination 
of [S* — x] + 1 independent exponential random variables with parameter 6, with coefficients 
X, X + 1, . . . , X + [S' — x] , if S* — X is not an integer. 

The n— th order moment of ^^(x) is given by (—1)" Mu^x)] ; calculating the first 
and second derivative, after some tedious computations, we obtain 

r^(x + 5-l) if ^-x = A;g1N 

^(^^(^))- j[5^(2x + [5-x]) if 5 -x^ IN 



k 

126*2 ' 



^^^j{12x2(A; + 1) + 12x(P - 1) + 3A;3 _ 2A;2 - 3A; + 2} if 5 - x = A; G IN 
mU^)) = I ^^{12x([S - x] + 2)(x + [5 - x]) + 

-[5-x](3[5-x]2 + 7[5-x] +2)} if^-x^lN 
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Remark 3.3 Notice that 



Ts{x) = inf{t >0:x + Nt>S} = mi{t > : S - x - Nt < 0} = ?{S - x) 

where r(y) is the first-passage time below zero of the process Y{t) := y — Nt. Moreover, it 
holds As{x) = Sts{x) — A{S — x), where A{y) is the area swept out by Y(t) till its first 
passage below zero. Thus, the moments of ts and the mean oi As can be obtained by the 
moments of r and E{A), which were calculated in [2]. In this way, one can avoid the heavy 
computations above. 



Example 5 (a Levy process) 

Let us consider the process X{t) = x + pt + Bt + Nt, where Nt is a homogeneous Poisson 
Process with intensity ^ > and let S* > 0, f3 + 9>S — x>0. We have ts{x) = mf{t > 
: X{t) > S} = inf{t > : {S - x) - /3t + Wt < Nt} i.e. the first hitting time of BM 



Wt with drift —/3 starting from S — x, and the Poisson process A''^ (see [TO] as regards the 
density of ts{x)) in a similar case). The condition (3 + 9 > S — x assures that ts{x) is finite 
with probability one (see Remark 12. 2p . Now, the infinitesimal generator of the process is 
Lf{x) = lf"{x) + (3f'{x) + f{x+ 1) — /(x) and the differential-difference equations involved 
to find the Laplace transforms of ts{x) and As{x), as well as those for the moments of ts{x) 
and As{x), cannot be solved explicitly; these quantities have to be found by a numerical 
procedure. 



Remark 3.4 Until now we have supposed that the starting point x < S" is given and fixed. 
One could introduce a randomness in the starting point, replacing X{0) = x with a random 
variable 77, having density g{x) whose support is the interval {—00, S). Thus, the quantities 
of interest become: Ts = infjt > : X(t) > S} and As = f^^ X{t)dt, while ts{x) and 
As{x) are the corresponding values conditional to t] = x. By using the found expressions for 
the moments of ts{x) and As{x), one can easily calculate the moments of ts and As- For 
instance, when X{t) is BM with drift fi, it follows by fl3.10p and fl3.14p that: 



E{A's) 



S — X 



g(x)dx = —E(S — rj) 



S — X {S — x) 



S — X 
2/i 



+ 



g{x)dx = -E{S - 77) + -E{{S - rif 



1\ 1 

S + X g(x)dx = — E 

2/i 



1^' 



(S-v)[S + v 



4 Conclusions and Final Remarks 

For S* > 0, we have considered a one-dimensional jump-diffusion process X{t) starting from 
X < S, that is, a diffusion to which jumps at Poisson-distributed instants are superimposed; 
then, we have studied the probability distribution of the (random) area As{x) = j^^^^^ X(t)dt 
swept out by X{t) till its first-passage time over the barrier S, i.e. the (random) time 
Ts{x) = inf {t > : X(t) > S\X{0) = x}. The analogous problem concerning the first passage 
of X{t) below zero was studied in [2], while results for Brownian motion with negative drift 
were obtained e.g. in pTO], [IB], pS], [SP, [23]. 
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In particular, we have shown that the Laplace transforms of As{x) and ts{x), their moments, 
as well as the probability distribution of the minimum of X{t) in [0,r5(a;)], are solutions 
to certain partial differential- difference equations (PDDEs) with outer conditions. Notice 
that, in the absence of jumps, these PDDEs with outer conditions become simply PDEs 
with boundary conditions. The quantities here investigated have interesting applications in 
Queueing Theory and in Economics (see the Introduction for a brief discussion). 

After considering theoretical results for diffusions X{t) with and without jumps, in the 
final part of the paper we reported some examples for which we have carried out explicit 
calculations. We remark that, in general it is not possible to solve explicitly the equations 
involved, in order to find closed formulae for the Laplace transform of the first-crossing area, 
and for its moments; moreover, even if one is able to find explicitly the Laplace transform, 
it is not always possible to invert it, to get the distribution of As{x). 

When the analytical solution is not available, due to the complexity of calculations, one 
can resort to numerical solution of the PDDEs involved; alternatively, one can carry out 
computer simulation of a large enough number of trajectories of the process X{t), in order 
to obtain statistical estimations of the quantities of interest. 
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